The source of elliptic flow and initial conditions for hydrodynamical calculations 
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^vj . A model for energy, pressure and flow velocity distributions at the beginning of relativistic heavy ion collisions 

is presented, which can be used as initial condition for hydrodynamical calculations. The results show that QGP 
forms a tilted disk, such that the direction of the largest pressure gradient stays in the reaction plane, but deviates 
from both the beam and the usual transverse flow directions. Such initial condition may lead to the creation of 

0\ \ "antiflow" or "third flow component" fl^. 
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►-j^ Introduction. Fluid dynamical models are widely used to describe ultra-relativistic heavy ion 

0^ \ collisions. Their advantage is that one can vary flexibly the Equation of State (EoS) of the matter and 
^^ ' test its consequences on the reaction dynamics and outcome. In energetic collisions of large heavy ions, 
Qs ■ especially if Quark-Gluon Plasma (QGP) is formed in the collision, one-fluid dynamics is a valid and 
<^ ' good description for the intermediate stages of the reaction. Here, interactions are strong and frequent, 
so that other models, (e.g. transport models, string models, etc., assuming binary collisions, with free 
propagation of constituents between collisions) have limited validity. On the other hand, the initial 
and final, Freeze-Out (FO), stages of the reaction are outside the domain of applicability of the fluid 
dynamical model. 

In conclusion, the realistic, and detailed description of an energetic heavy ion reaction requires a Multi 
P5 . Module Model, where the different stages of the reaction are each described with suitable theoretical 
approaches. It is important that these Modules are coupled to each other correctly: on the interface, 
L^ ■ which is a 3 dimensional hyper-surface in space-time with normal da^, all conservation laws should be 
'^ . satisfied (e.g. [T^^'^da^] = 0), and entropy should not decrease, [Sf^da^] > 0. These matching conditions 
P^ , were worked out and studied for the matching at FO in detail in refs. |j|]. 

The final FO stages of the reaction, after hadronization, can be described well with kinetic models 
where the matter is already dilute. 

The initial stages are more problematic. Frequently two or three fluid models are used to remedy the 
difficulties, and to model the process of QGP formation and thermalization. pi, 0, M Here, the problem is 
transferred to the determination of drag-, friction- and transfer- terms among the fiuid components, and 
a new problem is introduced with the (unjustified) use of EoS in each component in a nonequilibrated 
situations, where EoS does not exist. Strictly speaking this approach can only be justified for mixtures 
of noninteracting ideal gas components. Similarly, the use of transport theoretical approaches assuming 
dilute gases with binary interactions is questionable, as due to the extreme Lorentz contraction, in 
the C.M. frame, enormous particle and energy densities, with the immediate formation of perturbative 
vacuum should be handled. Even in most parton cascade models these initial stages of the dynamics are 
just assumed in form of some initial condition, with little justification behind. 

Our goal in the present work is to construct a model, based on the recent experiences gained in string 
Monte Carlo models and in parton cascades. One important conclusion of heavy ion research in the last 
decade is that standard 'hadronic' string models fail to describe heavy ion experiments. 
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All string models had to introduce new, energetic objects: string ropes ||5|, ^, quark clusters [^, fused 
strings R, in order to describe the abundant formation of massive particles like strange antibaryons. 
Based on this, we describe the initial moments of the reaction in the framework of classical (or coherent) 
Yang-Mills theory, following ref. |^ assuming larger field strength (string tension) than in ordinary 
hadron-hadron collisions. In addition we now satisfy all conservation laws exactly, while in ref. ||9[ infinite 
projectile energy was assumed, and so, overall energy and momentum conservation was irrelevant. We do 
not solve simultaneously the kinetic problem leading to parton equilibration, but assume that the arising 
friction is such that the heavy ion system will be an overdamped oscillator, i.e. yo-yoing of the two heavy 
ions will not occur. This assumption is based on recent string and parton cascade results. 

Formulation of model. Our basic idea is to generalize the model developed in |^], for collisions of 
two heavy ions and improve it by strictly satisfying conservation laws. First of all, we would create a 
grid in [x,y] plane (z - is the beam axes, [z,x] - is reaction plane). We will describe the nucleus-nucleus 
collision in terms of steak-by-streak collisions, corresponding to the same transverse coordinates, {xi,yj}. 
We assume that baryon recoil for both target and projectile arise from the acceleration of partons in 
an effective field F'^'^ , produced in the interaction. Of course, the physical picture behind this model 
should be based on chromoelectric flux tube or string models, but for our purpose we consider F^"^ as an 
effective abelian field. Phenomenological parameters describing this field must be fixed from comparison 
with experimental data. 

Let describe the streak-streak collision. 

d,T.Tr = T.Frn^,, (1) 

i i 

^mE< = 0' ^ = 1'2, (2) 

i 

nf is the baryon current of ith nucleus (we are working in the Center of Rapidity Frame (CRF), which 
is the same for all streaks. The concept of using target and projectile reference frames has no advantage 
any more). We will use the parameterization: 

nf = piuf , uf = (coshyj, sinhyj) . (3) 

T^^'^ is a energy-momentum flux tensor. It consists of five parts, corresponding to both nuclei and free 
field energy (also divided into two parts) and one defines the QGP perturbative vacuum. 

T,u = J2Tr + rpt = E [e^ ((l + Co) « - ch^") + n:i\ + Bg^'' , i = 1, 2 . (4) 

i i 

B ~ is the bag constant, the equation of state is Pi = c^ei, where e, and Pi are energy density and 
pressure of QGP. 

In complete analogy to electro-magnetic field 



Fr = d'^A^-d^A^=yi^ O^j , (5) 

Ui = d^A^i - d^Al , (6) 

/3 di^di^A^j 

U,i = -\n,uFr . (8) 



'"'f = "-i.T = 


-.p^{u'l±v^)=Pie^y^ , 


i i 

Tr+ T- — 





In our case the string tensions, cjj, will have the same absolute value a and opposite sign (in complete 
analogy to the usual string with two ends moving in opposite directions), and Uj will be constant in the 
space-time region after string creation and before string decay. 

To get the analytic solutions of the above equations, we use light cone variables 

(z,t) ^ (x+,x"), x^ = t±z. (9) 

Following [^, we insist that ei,yi, pi,A^ are functions of x~ only and 62,2/2,^25^2 depend on x+ only. 
In terms of light cone variables: 

(10) 

+ TF,^ , (11) 

where 

hi^ = (1 + cl)ei , /ii_ = (1 - cl)ei . (12) 

At the time of first touch of two streaks, t = 0, there is no string tension. We assume that strings are 
created, i.e. the sting tension achieves the value a at time t = to, corresponding to complete penetration 
of streaks through each other. 

Conservation laws — String rope creation. In light cone variables eq. (|2|) may be rewritten as 

(9_nj; + (9+n+ = . (15) 

So, we have a sum of two terms, depending on different independent variables, and the solution can be 
found in the following way. 

d^n^ = a, d+n^ = -a , .^g. 

n^ = ax^ + (ni)o, n^ = -ax^ + (722)0 ■ 

Since both n'^ and n^ are positive (and also more or less symmetric) we can conclude that for our case 
a = 0. 
Finally 

nr = pie-y' = poey° , n+ = ^26^^ = poe^° , (17) 

Pi = P^ey^^"' , P2 = Poe^«-^^ . (18) 

Let us come back to the energy- momentum tensor T^'^ . Based on eqs. (^, |^, ^) and taking into 
account the Lorentz gauge, d^A'^ — d'^Af = 0, we can find 

(T^,t T^,7\_( <rf-2a.{d^A^)+2a.{d^A^) ] nq^ 

[ T^; T^J )~\ a^ + 2a, {d^A^) + 2a, {d^A^) ) ' ^'^^ 

As mentioned before, after string creation, i.e. t > to, and before string decay we choose the string 
tensions in the form: 

a2 = -ai = a>0 . (20) 



To satisfy the above choice and the Lorentz gauge condition we take the vector potentials in the following 

form: 

At = 0, A^ = -2ax-, 

A+ = -2ax+, A^=0. ^^^^ 

In our calculations we used the parameterization: 

a = A(^] poVhf2, (22) 

where h, h are the lengths of initial streaks. The typical values of A are around 0.05 — 0.06. Notice, that 
there is only one free parameter in parameterization (p^). The typical values of a are 8 — 15 GeV/fm 
for £0 = 100 GeV/nud. 

The problem with eq. (|l|) is that we do not know what the really conserved quantities are. Using the 
definition of F'^'^ , eq. (y) , we can rewrite eq. (|l]) as 

d^T^^'^ = J2 Frni,>. = E (^^ (4"^*,m) - A'd^n,,^ - d^ (^fn,,^) + A^d^n,,^) . (23) 

i i 

The solution for n^ and 71-2^, eq. ([l7|), shows as that second and fourth terms vanish. So, we can define 
new energy- momentum tensor T^^'^, such that 

d^f'^'^ = , (24) 

t^- = J2 T^'"' + TpZt = E i^r - ^>t + r^>^o) + Bg^^ (25) 

i i 

Using the exact definition of J^ - eqs. (^) - we obtain 



rhuv _ i hi+e^y^ + 50-2 /ii_ + Aax'Tif \ ( h2+e'^y^ + cr^ + 2ax+n^ /12- - 2crx+n. 
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Now the new conserved quantities are 

Qo = Jf'^dV = Y:J f,''dV, (27) 

i ^ 

Q, = Jf'^'dV = J2j f,''dV. (28) 

Based on conservation of Qq, Q3 we can calculate rapidity, energy and baryon densities at the moment 
t = to, when the string with tension a is created. These new quantities are used as initial conditions for 
our differential eqs. (|^, |2|). 

^iM] = (Er\J: g^ I1+5I2 aey^ B h+h , . 

m ) \m) l+cl poeo(l+cg) ^h 8eo(l+cg)^' ^' Po£o(l+c§) 2/i ' ^ ' 

g2(to) \ ^ /£o\ _±_ _ tT^ ^2+5/i _ aeV" _ B h+h ,^^. 

m ) \m) l+cl /Ooeo(l+c§) ^h 8eo(l+cg)^' ^' Poeo{l+cl) 2h ' ^ ' 

Here the £i is energy per nucleon in CRF. Now the proper baryon density can be found, Piito) = Po— ffr) 

^. _ 1 _ £ifa) 
n 



^/W 



2 m 



For x^ > xq, where xq = 2to — |-2(0)| defines the string creation surface t = Iq, for nucleon or cell 
element in the position z = z{0) at the time t = 0, we should solve eqs. (p^), with boundary conditions 

Hi {x~ = Xq) = pQC^y^ n2 {x~^ = Xq) = /9oe^^° 

hi+ix" = Xq) = ei(to)(l + Co) h2+ix+ = Xq) = ei(to)(l + Co) 

yi(x~ = Xq) = yi(to) y2{x^ = Xq) = ^2(^0) 

ai{x~ = Xq) = -a Cr2(x+ = Xq) = a , 

where ej(to) = fTiPiito) ~ energy density in the rest frame of ith nuclei. 
Let us present the complete analytical solution in the following form 



,(-r+^2y. 



di fdi 



hi \hi 



ii ^(-r+'2j,,(to) 



X — Xq 



OLaA 



(31) 



(32) 



h+ = e^-y^"^y^ei{tQ){l + cl)e 



,2v-(-)'+'23/,(to) 



X — Xq 



where x^ = x~ , x^ = x+, i,j = 1,2 , i j^ j, and using the notations 

bi = aaj + 2(T/Ooe^° , 
d, = Q - 2cT/>oe^Oe(-)'^'2y.(io) ^ 
ei(to)(l + cg) 



Ti 



,(-)'+i2j/,(io)-- 



ai = ci + 2apoe^o - 2apoe^Oe22^i(*o) , as = C2 + 2a/9oe^» - 2apoe^°e-2?/2(*o) 

Ci = q((1 + Co)ei(to) - eo)/2to . 
Then the trajectories of nucleons (or cell elements) for both nuclei are given by: 



x+{x-) = z{0) + J^^ dx e2?/i(^) 






^(0)-f7i2; -XQ 



) + (fr 
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n 
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2apQeyO 
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-1 


2crpoeyo 



■(x+) = -z(0) + /,f dx e-2s^2W 



(X^ — Xo 



)+( 



d2 



-z{^)-^\^ ^ -^u; -r \^ ft^ -r o - ■; '^2apoe''0 



_^g-2j/2(to)] -j-^ ^ ""1 



2 / 



J"POe 



for nucleon or cell element in the position z = z(0) at the time t = 0. 



(33) 
(34) 

(35) 
(36) 

(37) 

(38) 
(39) 



(40) 



(41) 



Recreation of the matter. As we may see from the trajectories, eqs. ( ^ 41), nucleons (or cell 
domains) will keep going in the initial direction up to the time t = ti^tum, then they will turn and go 
backwards until the two streaks again penetrate through each other and new oscillation will start. Such 
a motion is analogous to the "Yo-Yo" motion in the string models. Of course, it is difficult to believe 
that such a process would really happen in heavy ion collisions, because of string decays, string-string 
interactions, interaction between streaks and other reasons, which are quite difficult to take into account. 
To be realistic we should stop the motion described by eqs. (4^, 41) at some moment before the projectile 
and target cross again. 

We assume that the final result of collisions of two streaks after stopping the string's expansion and 
after its decay, is one streak with length Alf, homogeneous energy density distribution, e/, and baryon 
charge distribution, pj, moving like one object with rapidity yj. We assume that this is due to string- 
string interactions and string decays. As it was mentioned above the typical values of the string tension, 
a, are of the order of 10 GeV/ fm, and these may be treated as several parallel strings. The string-string 
interaction will produce a kind of "string rope" between our two streaks, which is responsible for final 
energy density and baryon charge homogeneous distributions. Now it is worth to mention that decay of 
our "string rope" does not allow charges to remain at the ends of the final streak, as it would be if we 
assume full transparency. 

The homogeneous distributions are the simplest assumptions, which may be modified based on ex- 
perimental data. Its advantage is a simple expression for ej, pf, yj. 

The final energy density and rapidity, e/ and yf, may be determined from conservation laws. 



{M\l + eg) + 2cg) + ^{M^l + eg) + 2c2)2 + 4c^(M2 - 1) 

COSh^'w/- = ^ rj- — ^ ^ , (42) 

^^ 2(l + cg)(M2-l) ' ^ ^ 

where we neglected B A If next to Qq and introduced the notation M = {h + h)/{h — h), 

■7^-BAlf 

er = ^^^ i (43) 

^ {{I + cl) cosh' yf- 4) Alf ■ ^ ^ 

The typical trajectory of the streak ends is presented in Fig. [l|. From Iq they move according to eqs. 
|, 41) until they reach the rapidity yi = yf- Later the final string starts to move like one object with 



rapidity yf. 

The turning points can be found from the condition: 

yi = yf , (44) 

which gives for ith. nucleus (xi = x~ , X2 = x"*") 



Xi turn — Xq + Tj 



1 






(45) 



Initial conditions for hydro dynamical calculations. In this section we present the results of 
our calculations. We are interested in the shape of QGP formed, when string expansions stop and their 
matter is locally equilibrated. This will be the initial state for further hydrodynamical calculations. We 
may see in Figs. |2|, that QGP forms a tilted disk for b ^ 0. So, the direction of fastest expansion, the 
same as largest pressure gradient, will be in the reaction plane, but will deviate from both the beam axis 
and the usual transverse flow direction. So, the new flow component, called "antifiow" or "third flow 
component" , will appear in addition to the usual transverse flow component in the reaction plane. With 
increasing beam energy the usual transverse flow is getting weaker, while this new flow component is 
strengthened. The mutual effect of the usual directed transverse flow and this new "antifiow" or "third 
flow component" lead to an enhanced emission in the reaction plane. This was actually observed and 
widely studies earlier and referred to as "elliptic flow". 
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Figure 1: The typical trajectory of the ends of two initial streaks, corresponding to numbers of nucleons, 
ni and n2- Stars denote the stopping and turning points, where yi = yf. From to to the turning points 
streak ends keep going in their initial direction according to eqs. ( ^0| , 41). Later the final streak starts 
to move like one object with rapidity, yf ( |4^ in CRF. 



Conclusions. Based on earlier Coherent Yang-Mills field theoretical models, and introducing effective 
parameters based on Monte-Carlo string cascade and parton cascade model results, a simplified model is 
introduced to describe the pre fluid dynamical stages of heavy ion collisions at the highest SPS energies 
and above. The model predicts limited transparency for massive heavy ions. 

Contrary to earlier expectations, — based on standard string tensions of 1 GeV/fm which lead to the 
Bjorken model type of initial state, — effective string tensions are introduced for collisions of massive 
heavy ions, as a consequence of collective effects related to QGP formation. These collective effects 
in central and semi central collisions lead to an effective string tension of the order of 10 GeV/fm and 
consequently cause much less transparency than earlier estimates. The resulting initial locally equilibrated 
state of matter in semi central collisions takes a rather unusual form, which can be then identified by 
the asymmetry of the caused collective flow. Our prediction is that this special initial state may be the 
cause of the recently predicted "antiflow" or "third flow component". 

Detailed fluid dynamical calculations as well as flow experiments at semi central impact parameters 
for massive heavy ions are needed at SPS and RHIC energies to connect the predicted special initial state 
with observables. 
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Figure 2: The Au+Au collisions, eo = 100 GeV/nud, b = 0.5{Ri + R2), A = 0.055 (parameter A 
introduced in (|2^)), y = (ZX plane through the centers of nuclei). We would like to notice that final 
shape of QGP volume is a tilted disk ~ 45^, and the direction of the fastest expansion will deviate from 
both the beam axis and the usual transverse flow direction, and might be a reason for the third flow 
component, as argued in p^]. 



K 

[5 

[6: 

[7; 

[9 
[10: 



L.P. Csernai et al., Phys. Rev. C26 (1982) 149. 

J. Brachmann, S. Soff, A. Dumitru, H. Stoker, J. A. Maruhn, W. Greiner, D.H. Rischke, L. Bravina, 
Phys. Rev. C61 (2000) 024909. 

T.S. Biro, H.B. Nielsen, J. Knoll, Nucl. Phys. B245 (1984) 449. 

K. Wegner, J. Aichelin, Phys. Rev. Lett. 76 (1996) 1027. 

N.S. Amelin, M.A. Braun, C. Pajares, Phys. Lett. B306 (1993) 312, Z. Phys. C63 (1994) 507. 

H. Sorge, Phys. Rev. C52 (1995) 3291. 

M. Gyulassy, L.P. Csernai, Nucl. Phys. A460 (1986) 723. 

L.P. Csernai, D. Rohrich, |nucl-th/990803^ , Phys. Lett. B458 (1999) 454. 



